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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


RESEARCH MEMORANDUM 


EFFECT OF SOME SECTION MODIFICATIONS AND PROTUBERANCES 
ON THE ZERO-LIFT DRAG OF DELTA WINGS AT TRANSONIC 
AND SUPERSONIC SPEEDS 


By Cari A. Sandahl and William E. Stoney 
SUMMARY 


Free-flight rocket model tests have been made at Mach numbers 
extending from approximately 0.8 to 1.9 to determine the effect of 
leading-edge radius апа airfoil aftersection thickening on the zero-lift 
drag of a 60° delta wing. In addition, one configuration employed simu- 
lated control~actuator housings. For the subsonic leading-edge condi- 
tion of the tests, increasing the leading-edge radius had no measurable 
effect on the drag coefficient. A small but finite thickening of the 
airfoil aftersection which did not increase the trailing-edge thickness 
resulted in a slight increase in drag coefficient. Further thickening 
of the aftersection resulted in substantial increases in drag coefficient. 
At supersonic speeds, the drag coefficient increased nearly linearly with 
those values of trailing-edge thickness tested. Increasing the section 
thickness ratio near the wing tip while maintaining straight-line 
constant-percent-chord wing~surface elements resulted іп no measurable 
drag increment. Linear theory generally predicted absolute values of 
wing pressure drag coefficients larger than those derived from the experi- 
ments, but predicted effects of section changes accurately. 


INTRODUCTION 


Thin deita wings аге of considerable current interest in the design 
of transonic airplanes. Іп such designs it often appears desirable from 
structural and volume considerations to increase the wing section thick- 
ness fore апа aft of the position of maximum thickness while maintaining 
the maximum thickness constant. It may also be desirable to increase 
the leading-edge radius in order to reduce the drag at lift. Generally, 
all the foregoing section changes may be expected to increase the zero- 
lift drag. It was the primary purpose of the investigation reported 
herein to assess the zero-lift drag penalties arising from such section 
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modifications for а 609 delta wing. The basic airfoil section was of the 
МАСА 000Х-Х5 family of reference 1. The leading-edge radius was varied 
from zero to three times that normally used with the basic section. The 
airfoil aftersection was varied systematically from that of the basic 
section to'one having parallel sides extending back from the maximum- 
thickness position. For one test, the basic wing was equipped with 
rather large simulated control-actuator housings. 


The Mach number range of the tests was from approximately 0.8 to 1.9. 
The corresponding Reynolds numbers varied from 6.5 x 10° to 27.4 x 106 
based on the wing mean aerodynamic chord. 


Free-flight rocket-propelled test vehicles were employed. The 
models were launched at the Langley Pilotiess Aircraft Research Station 
at Wallops Island, Va. 


SYMBOLS 
Съ configuration drag coefficient based оп 5, 
CD wing-plus-interference drag coefficient based on Sy 
Орр wing pressure drag coefficient based оп Sy 
бү wing area obtained by extending leading and trailing edges to 


center line (equal to 7.57 sq ft for 811 models except 
model 10 for which бу equals 4.55 sq ft) 


W model weight, sustainer motor fuel expended 

ү velocity tangent to flight path or volume of wing external to 
fuselage - 

а, linear acceleration tangent to flight path 

& acceleration due to gravity E 

e mass density of air Е 

y flight-path angie measured from horizontal _ | == 

М Mach number 

R Reynolds number based on С 
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Атр leading-edge sweep 
Атр trailing-edge sweep 
b actual wing span 
b! wing span measured to theoretical wing tips 
с local wing chord 
Ср wing chord over which section sides аге parallel 
с wing mean aerodynamic chord 
t Section maximum thickness 
h section trailing-edge thickness 
A configuration cross-sectional are& 
г radius of equivalent body of revolution 
x distance from fuselage nose 
1 fuselage length 
MODELS 


The gener&l arrangement of the test configuration is shown in fig- 
ures 1 and 2. The fuselage ordinates are presented in table Т. A list 
ОР configurations tested is in table II. 


The fuselages were constructed mainly of laminated wood. The wings 
had а balsa-mahogany core to which 0.020-inch-thick duralumin plates were 
bonded to form the outer surfaces. The wing tips and trailing edges were 
stiffened by means of a chord-plane insert of 0.064-inch-thick duralumin 
The vertical tails were of duralumin. 


The geometry of the wings used is shown in figures 3 to 6. The air- 
Poil sections are defined in reference 1. The wing of model 11 was a 
60° delta wing with ап NACA 65А005 airfoil section. The notched wing 
of model 10, figure 5, was obtained by removing the flat region of the 
wing of model 6, figure 4, and shearing the airfoil aftersections for- 
ward. Both wings had the same absolute maximum thickness. 


4 "OONFIDENBEAT, NACA RM L52L2Àla 


The wing of model 12 was equipped with simulated control-actuator 
housings as shown in figure 6. The housings were scaled from those on 
the top side of the wing of a current airplane. For the test, duplicate ` a 
housings were employed on both sides of the wing to ensure zero-lift— 
flight and а measurable drag increment. The drag-coefficient increment 
would be approximately doubled by this procedure. The housings were of 
maximum breadth and thickness at, and were all identical aft of, the 
hypothetical control hinge line which was located at 90 percent of the 
theoretical root chord and was parallel to the wing trailing edge. In 
side view, the housings are defined by parabolic arcs with vertices 
above and below the chord plane at the hinge line, tangent to the wing 
surface at the front end, and intercepting the trailing edge. In plan 
view, the housings are defined by similar parabolic arcs. А photograph 
of this model is shown in figure 7. 


The longitudinal distribution of cross-sectional_area for all models 
is shown in figures 8 to 10. | 


TEST METHOD 
Propulsion 


The models were accelerated to supersonic speeds by a two-stage Е 
rocket propulsion system. The first stage was either а 5-inch НРАС 
rocket motor or а 6-inch ABL Deacon rocket motor yielding maximum Mach 
numbers of approximately 1.5 and 1.9, respectively. Тип all cases the 
second stage (sustainer) was а 3.25-inch Mk. 7 rocket motor. А photo- 
graph of the model-booster-launcher arrangement preparatory to launching 
is shown in figure 11. = = 


Instrumentation 


The instrumentation, none of which was contained within the models, 
included a CW Doppler velocimeter radar, a modified SCR 584 position 
radar, and radiosonde equipment. The Doppler velocimeter produced a _ 
time history of the radial component of the model velocity relative to 
the launching site. The position radar produced a time history of model 
position in space. The radiosonde units, which were launched immediately 
after the test flights, were tracked during their ascent with the space 
position radar permitting the determination of the speed and direction 
of winds aloft in addition to the usual pressure and temperature i 2 == 
measurements. 5 Do 


МАСА RM L53L2ha дами? AT. 5 


Рафа, Reduction 


The drag coefficients of the models were obt&ined during coasting 
flight by evaluating the expression 


-2W(a + в sin у) 


Бру Sy 


Со = 


In the evaluation, the position-radar record was reduced to a plot of 
altitude against horizontal range. From this flight-path plot, trigono- 
metric corrections were obtained which were applied to the CW Doppler 
data to convert radial velocity to tangential velocity along the flight 
path. The component of the wind velocity parallel to the flight-path 

was then added to obtain airspeed. The acceleration along the flight 
path was obtained by numerical differentiation of the tangential velocity. 
Mach number and Reynolds number were also obtained from the radar and 
radiosonde data. 


The errors are estimated to be within the following limits: 


Ср e е о э ө е е а е » ее в ө Қ e è> ее е е е е 9 ооо в ө £0 .0007 
М “ее c£? чов е ае ө е © o © ө ùo е $9 ә ù ее а ò е э в е +0 .005 


The errors in wing-plus-interference drag coefficients obtained by sub- 
tracting fuselage drag from wing-fuselage drag may be somewhat larger. 
A typical set of test results is shown in figure 12 to illustrate the 
continuity and scatter of the data. 


RESULTS AND DISCUSSION 


The variation of Reynolds number with. Mach number for the tests is 
shown in figure 13. The results for each model are presented in figure 14 
wherein are plotted the total-drag coefficients and the wing-plus- 
interference drag coefficients obtained by subtracting body-drag coef- 
ficients from the total-drag coefficients. The experimentally determined 
body-drag coefficients, based on the two values of total wing area 
employed in the tests, are shown in figure 15. The curve labeled "4 fins" 
was obtained from & free-flight test of & four-fin wingless body identi- 
cal to that used with the winged models of the tests. The decrement in 
drag resulting from decreasing the number of fins from four to two was 
obtained from experimental results in reference 2 obtained with geometri- 
cally similar fins mounted on а cylindrical fuselage. The drag coeffi- 
cients for the body with the two fins so obtained were subtracted from 
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the drag coefficients for the winged two-fin models to obtain wing-plus- 
interference drag coefficients. It should be noted that this subtrac- 

tive process renders the wing-plus-interference drag coefficients sub- un 2-2 
Ject to additional errors in the regions where the drag coefficients қ 
change rapidly with Mach number (as near M = 1.0) because of the сош- = 
bined errors in Mach number and drag coefficient. | 


Effect of Leading-edge rađius.- The effect of leading-edge radius - 
is shown in figure 16 to be small and within the experimental accuracy 
for the range of radii tested. It should be noted that the leading edge | 
was always subsonic. = 


Effect of flat midsection.- The effect of flat midsections is shown 
in figure 17. No measurable increases in drag coefficient over that of 
the basic section were obtained for с/с = 0.2. The larger values of 


сг/с resulted in substantial reductions іп critical Mach number and 
increases in supersonic drag. For ср/с = О.Т, the trailing-edge thick- 


ness becomes equal to the maximum section thickness. Pressures acting 
on this base cause the large drag increment obtained for this configura- 
tion over the entire Mach number range investigated. Р 


In figure 18 is shown а comparison of wing pressure drag coeffi- 
cients calculated by the method described in the appendix with those е 
derived from the даља Рог the models of figure 17. The wing pressure 
drag coefficients were derived from the experimental results by sub- 
tracting the wing-friction-drag coefficients from the wing-plus- 
interference drag coefficients at the Mach numbers noted in figure 18. 
The wing-friction-drag coefficients were obtained by adjusting the 
subsonic-wing-plus-interference drag coefficient of the basic configu- 
ration (model 3) to the appropriate Mach and Reynolds numbers using the 
theoretical variations of turbulent-skin-friction coefficient with Mach 
number and Reynolds number given in reference 5. In figure 18(а) it 
will be noted that the theoretical drag coefficients are generally larger 
than those experimentally derived. In figure 18(b) the theoretical and 
the derived increments in pressure drag coefficients due to the flat 
midsections are seen to be in good agreement. It should be noted that en 
the theory increasingly overestimates the pressure drag coefficient as 
the Mach lines approach the leading edge. _ 
In figure 19 is plotted ihe ratio of wing-plus-interference drag 
coefficient of the flat midsection wings to that of the basic section m 
as a function of the relative wing volume. The curves were faired by 
giving less regard to those test points at а relative volume of 1.11 _ 
because they арреаг questionable on the basis of theory noted previously. 
The results indicate that some increases in wing volume can be obtained я 
with small drag penalty. Further increases in volume result in large 
increases in drag. 
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Effect of trailing-edge thickness .- The effect of trailing-edge 


thickness is shown in figure 20 wherein are plotted the total-drag coef- 
ficients for configurations having trailing-edge thicknesses equal to О, 
30, and 100 percent of the wing maximum thickness. At supersonic speeds, 
the drag coefficients increase nearly linearly, within the experimental 
accuracy, with those values of trailing-edge thickness tested. Also 
shown in figure 20 are drag-coefficient curves obtained, at supersonic 
speeds, by adding the increments in wing drag coefficient due to trailing- 
edge thickness, calculated as described in the appendix, to the measured 
drag coefficients of model 5 (h/t - 0). In the Mach number range from 0.8 
to 1.0, the calculated values were obtained by adding the base drag coef- 
ficients, calculated from the base pressure measurements of reference |, 
to the drag coefficients of model 5. Excellent agreement is shown in 
this comparison at supersonic speeds. At subsonic speeds, the calculated 
values are substantially higher than the experimental values, indicating 
@ reduction in wing drag exclusive of that contributed by the base with 
increasing h/t. In figure 21 is shown a comparison of the wing pres- 
sure drag coefficients calculated by the method of the appendix with 
those derived from the experiments in the manner described in the dis- 
cussion of figure 18. The calculated values are generally higher than 
those experimentally derived; the calculations do, however, indicate 
accurately the trends. 


Effect of trailing-edge sweep.- In figure 22 is shown the effect 


on the wing-plus-interference drag coefficient of removing the flat 
region of wing of model 6 (ce/c = 0.4) and shearing the aftersections 
forward. The thickness ratio of the resulting wing varied from 0.05 at 
the theoretical root to 0.10 at the 90-percent semispan station. This 
wing was considered to be structurally feasible while offering the ров- 
sibility of low drag at low supersonic speeds. The drag coefficients 
are based on the total wing areas of the individual models. Substantial 
drag reductions were indeed obtained at Mach numbers less than approxi- 
mately 1.2; at higher Mach numbers, the sheared wing had higher drag. 

As noted previously, the abrupt changes in CD. near М = 1.0 are 


believed to result from the subtractive process by which this quantity 
was obtained. 


Effect of model scale.- The effect of model scale is illustrated 
in figure 25 wherein the results for model 11 are compared with those 
for а geometrically similar twice scale model (model 4, from ref. 2). 
The differences are within the experimental errors and, particularly 
for total configuration drag coefficients at supersonic speeds, are in 
the direction and of the magnitude calculable on the basis of changes 
in skin-friction coefficient. 


Effect of spanwise variation of section thickness ratio.- The effect 


of spanwise variation in thickness ratio may be noted in figure 14 by 
comparing the results for model 3, in which the wing section thickness 
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ratio varied from 0.03 at the root to 0.06 at the 90-percent-semispan 
station, with those for model 11 in which the section thickness ratio 
was constant at 0.03. Calculations by the method of the appendix showed 
that no measurable differences in wing drag could be attributed to the 
different wing sections employed in-the two models. For the range of 
Mach numbers for which the foregoing comparison can be made, no measur- 
&ble differences in drag exist. 7 “ 


Effect of simulated control-actuator housings.- The effect of the 


control-actuator housings on the drag is shown in figure 24 wherein are 
compared results for the basic wing, the wing with housings, ала а wing 
with a flat airfoil midsection with approximately the same volume ав 
that of wing and housing. At low supersonic speeds, the increment in 
drag coefficient is substantial and approximately as large as that for Е 
the wing having с/с = 0.4. It should be noted, as pointed out in the 

"Models" section, that this increment is approximately twice that which 

would be obtained for & more realistic situation with the housings on 

only one side of the wing. Also shown in figure 24 is the estimatéd  . ©. 
drag coefficient for model 12 obtained by adding to the drag coefficient 
of model 5 an increment estimated for the housings by (1) calculating 
the longitudinal distribution of cross-sectional area external to the 
basic wing produced by the upper апа lower housings &t each spanwise 
station &t which housings were attached, (2) obtaining the equivalent 
bodies of revolution (6 in number), and (5) estimating (from ref. 5) " 
and summing the drag of the bodies. The estimate so obtained is low at 

low supersonic speeds because of interference effects between the housings 

and improves with increasing Масћ number as these interferences become 

smaller. 


CONCLUSIONS 


The present investigation made to determine the effect of some modi- 
fications and protuberances on the zero-lift drag of delta wings at 
transonic and supersonic speeds indicated the following: 


1. For the subsonic leading-edge condition of the tests, increasing 
the leading-edge radius had no effect on the drag coefficient. 


2. A small but finite thickening of the airfoil aftersection which 
did not increase the trailing-edge thickness resulted in no measurable 
increases in the drag coefficient. Additional thickening of the after- E 
section resulted in substantial increases in drag coefficient. 


2. At supersonic speeds, the ‘drag coefficient increased nearly 
linearly with those values of trailing-edge thickness tested. 


didam, 
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4. Increasing the wing section thickness ratio at the 90-percent- 
semispan station from 0.03 to 0.06 while maintaining straight-line 
constant-percent-chord wing-surface elements and a root thickness ratio 
of 0.05 resulted in no measurable drag increment. 


5. Multislope linear-theory methods for calculating wing pressure 
drag coefficients combined with experimental base pressure values gen- 
erally yielded &bsolute values which were larger than those derived, from 
the experiments, but accurately predicted effects of section changes. 


Langley Aeronautical Laboratory, 
National Advisory Committee for Aeronautics, 
Langley Field, Va., December 8, 1953. 
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APPENDIX 
THEORETICAL CALCULATIONS 


Ihe theoretical calculations of the wing wave drag were made by the 
method of reference 6. In this method the airfoil section is approxi- 
mated by а many«sided polygon and the drag coefficients are obtained by 
superimposing the known solutions for single triangular wings of con- 
stant slope. This is identical to the procedure followed in reference. 7 
in the calculation of the drag of two- and three-slope wings. 


The number of sides necessary depends somewhat on the shape of the 
airfoil section; however, it was shown in reference 6 for a circular- 
arc wing that the calculated drag was essentially constant above a divi- 
sion into 10 slopes. One of the present airfoils was_calculated using 
lO-sided and 20-sided figures and the results were identical for Mach 
numbers where the leading edge was subsonic. For cases where the Mach 
line lies оп or behind the leading edge, the дгар calculated by this 
method for any round-nose airfoil is incorrect since it approaches infin- 
ity as the number of sides is increased. This is due to the stagnation 
condition at the leading edge for which the linear theory is inapplicable. 
Ihe present calculations were extended to the supersonic leading-edge 
cases since it is the change in drag with change in shape after the maxi- 
mum thickness that is of interest and this should be reasonably independ- 
ent of nose shape. Although none of the experimental results include 
the supersonic leading-edge condition, the calculated results are included 
as a matter of interest. Е 


The results of these calculations are shown in figures 25 апа 26 
for wings having a uniform thickness ratio of 0.03. The results in fig- 
ure 25 show the effect of the extent of flat airfoil midsection. The 
corresponding test models are noted. The results in figure 26 show the 
effect of decreasing the average airfoil aftersection slopes from those 
of the basic NACA 0003-63 section. The values of h/t resulting from 
this procedure are noted in the figure. For those cases in which a 
finite thickness at the trailing edge existed, a base drag allowance 
was included utilizing the base pressure measurements from reference 8. 
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TABLE ТТ 


TESI CONFIGURATIONS 


NACA airfoil section 
А Arr = 
Кобе eee or [ang ее ЫЈ) | Theoretical root __|0-9 theoretical semtspan | 
ratio |deg | deg ег/ / Pt Theoretical root 0.9 theoretical semispan 
О to 0.3c | 0.3c to 1.06 | О to 0.5с | 0.5е to 1.0с 
2.06 | 60 


0003-03 0006-03 
0003-33 0006-33 
0003-63 0006-63 








0003-93 0006-93 


0005-65 See fig. | 0006-65 See fig. 
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See fig. 4 
See flg. 4 
See fig. 4 
See fig. № 


oee fig. 5 


654005 


0005-65 


0006-65 
0006-65 
0006-65 
0006-65 


бее fig. 
See fig. 
pee fig. 


See fig. 


See fig. 5 
654005 
0006-65 





ет 








6.50 Q487 Рат plate, 
beveled КАТИЋ and 


trawling edges 


Figure ].- General arrangement. Dimensions are in inches. 
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Figure 2.- Photographs of typical model. 
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Figure 5.- Geometry of wings used in investigation of leading-edge effects, 
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_ Figure 4.- Geometry of wings used in investigation of trailing-edge effects. 
Numbers denote models. 
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Figure 5.- Geometry of notched delta wing. Model 10. 
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Figure 6.- Geometry of wing of model 12. Dimensions are in inches. 


те 


oono 





NACA ™ L55L2ha 
20 


Figure Т.- Model 12. 
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Figure 8.- Longitudinal distributions of cross-sectional area. Numbers 
denote models. 
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Figure 9.- Longitudinal distributions of cross-sectional area. Numbers 
denote models. | 
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Figure 10.- Longitudinal distribution of cross-sectional area. Model 12. 
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Figure 11.- Model~booster-launcher assembly. 
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Figure 12.- Typical set of test results. 
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Figure 14.- Test results. 
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Figure Щ.- Continued. 
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(h) Model 8. 


Figure 14.- Continued. 


50 | Енреа МАСА ЕМ L535L2ha 


OL : : TE Hi F 
FHH ШЕ НН 
Е 


ЕА НЕ 
ЖЫНЫ 1 : | үүн. 





S /O /2 #4 LO LE ZO 
^7 


(1) Model 9. 
ее НЕНА НЕНА E | ЕНЕ 


- + TER LLELTTT A HERE E HF 
ПНЕ 


22 Ни ЕЕ ее еее рев реве еве 


CD Ol 





(k) Model 11. 


Figure 14.- Continued. 
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(1) Model 12. 


Simulated control-actuator housings. 


Figure 14.- Concluded. 
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Figure 15.- Body dreg. 
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Figure 16.- Effect of leading-edge radius. 
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Figure 17.- Effect of flat extent. 





NACA ЕМ L53L2la dE 55 


‚989 





об 
Ў 
5 
5 A 
N 
NI 
J 
А 008 
Х\ 
„OOZ 
2, 
(a) Total drag coefficients. 
a 
3 
S 
S 
55 
5 
Q 
oe 
N 





ACp, (From fests) 


(b) Incremental drag coefficients measured from those 
of the unmodified section (model 3). 


Figure 18.- Comparison of experimentally derived and calculated pressure 
drag coefficients for series of wings in which extent of flat mid- 
section is varied. Numbers denote models. 


POOR ыы. 
е с=т= 


56 eU ES МАСА RM L55L2ha 





2 >” La 42 IL Да 
Volume (With Hat) / Volume (Basic) 


Figure 19.- Variation of relative wing drag with relative wing volume 
for series of wings in which extent of flat midsection is varied. 
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Figure 20.- Effect of trailing-edge thickness. 
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Figure 21.- Comparison of experimentally derived and calculated pressure 
drag coefficients for series of wings in which trailing-edge thickness 


is varied. Numbers denote models. 





Figure 22.- Effect of shearing afterpart of wing forward. 
same absolute thickness. 
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Figure 25.- Effect of model scale on drag. 
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Figure 23, Effect of simulated control-actuator housings on drag. 
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Figure 25.- Calculated pressure drag coefficients for series of wings in 
which extent of flat midsection is varied. t/c = 0.03. 
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Figure 26.- Calculated pressure drag coefficients for series of wings in 
which aftersection slopes аге varied. %/с = 0.05. 
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